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We study the effect of superconducting fluctuations on the tunnel current-voltage characteristics 
of disordered superconducting films placed in a perpendicular magnetic field, H , in the whole H-T 
phase diagram outside the superconducting region. This tunnel-current is experimentally accessible 
by STM measurements. In the domain of temperatures T > T c q and relatively weak fields H <C 
H C 2(0) we reproduce existing results for the zero- voltage tunneling conductance, but also discover 
an important nonlinear contribution, which appears due to dynamic fluctuation modes and results 
in the formation of a strong zero-bias anomaly (ZBA) on the scale eV ~ £;b(T — T c q). At large 
voltages (eV ~ feTco) these modes, together with the contribution from static fluctuations, form 
a pseudogap maximum. At low temperatures, with magnetic field values near H C 2(0), fluctuations 
acquire quantum character and the general picture of the voltage dependent tunneling conductance 
resembles that one close to T c o, where the role of temperature and magnetic field are exchanged. In 
particular, a gap-like structure appears with maximum at eVi na x ~ Abcs and a sharp ZBA on the 
scale eV ~ Abcs(H / H C 2(0) — 1). The complete expression for the tunneling current at arbitrary 
fields and temperatures can be evaluated only numerically, which is presented in detail. 



PACS numbers: 74.40.-n 

I. INTRODUCTION 

According to the microscopic BCS thcoryS, the super- 
conducting state is characterized by a gap in the normal 
excitation spectrum, centered at the Fermi level, Ep, 
which vanishes along the transition line Hr*{T). How- 
ever, it was predicted, as early as in 1970a, that even 
above the superconducting critical temperature, T c q in 
zero magnetic field, in its normal phase, thermal fluc- 
tuations result in a noticeable suppression of the den- 
sity of states, v(E), of the superconductor in an en- 
ergy range around the Fermi level. rMore specifically, 
in the case of a disordered thin filmELthe characteris- 
tic energy range of this suppression isa Eq ~ T — T c o 
and its zero energy depth is given by (E = 0,T) = 

-0.1Gi (2) i/„[T c0 /(T-T c0 )] 2 , where Gi (2) = l.3/(p 2 F ld) 
is the Ginzburg-Levanyuk number characterizing the 
strength of fluctuations in the film, and v n is the den- 
sity of the states of a normal metal at the Fermi level. 

With the discovery of high temperature superconduc- 
tivity, the investigation of the suppression of the density 
of states and the domain of parameters, where it is ob- 
served in cuprate superconductors developed into a major 
research area. It became known as the pseudogap state 
and is believed to be a key mechanism for high tempera- 
ture superconductivity itself]. A recent experimental re- 
port of the pseudogap in a conventional 2D superconduc- 
tor, a disordered ultrathin titanium nitride (TiN) filmQ, 
added another piece of the pseudogap puzzle. Com- 
bined scanning tunneling microscope (STM) and trans- 
port measurements of the I-V characteristics allowed to 
carefully compare these experimental findings with the- 



oretical predictions of the behavior of the contribution 
to transport conductivity due to superconducting fluc- 
tuations (SF) [0. As a result, the observed pseudogap, 
existing over a wide range of temperatures above T c q, was 
attributed to two-dimensional (2D) SF, which are favored 
by the proximity to the superconductor-insulator transi- 
tions. This recent finding paralleled an earlier observa- 
tion of the pseudogap due to SF in experiments on Al-I-Al 
tunnel junctions!! that have been observed soon after its 
prediction and poses a call to re- visit the pseudogap state 
in high-T c cuprates experimentally and carefully inspect 
the role of superconducting fluctuations. 

A major experimental tool for determining the density 
of states is by measurements of the differential tunnel 
conductivity 

<rtun (V) = dl (V) /dV 

~/:(-^)"< e +- ><*■ <» 

Considering that the success in revealing the nature of 
the pseudogap in superconducting filmsEl was based on 
the ability to identify the effect of SF on the electronic 
transport, one immediately recognizes the quest for un- 
covering the role of fluctuations in tunneling properties. 
The knowledge of the behavior of the fluctuation contri- 
bution to er tU n would have offered an irreplaceable tech- 
nique for identifying fluctuation effects directly from the 
experimental data. Our work addresses this challenging 
task and develops a theory for the fluctuation contribu- 
tion to the tunneling conductance. 
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II. MODEL 

We study the effect of SF on the tunneling current 
I (V) between a normal metal electrode and a disordered 
two-dimensional superconducting film placed in perpen- 
dicular magnetic field throughout the whole phase dia- 
gram above the H C 2(T) line. Describing this system by 
means of a tunnel Hamiltonian, the tunnel-current can be 
expressed in terms of the correlator K (w„) of the elec- 
tron Green's functions of the corresponding electrodes, 
which is analytically continued from Matsubara frequen- 
cies ui v = 27r7V, v = 0,1,2,... to the upncr half-plane of 
complex frequencies u„ — > —iuj = —ieV}i: 



<j(V) = ~e!mK R {eV). 



(2) 
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FIG. 1. (Color online) Diagram contributing to the fluctu- 
ation tunnel-current in first orders of barrier transparency 
and strength of fluctuations Gi(2) . Solid lines correspond to 
single-electron Green's functions of the electrodes, the wavy 
line represents the fluctuation propagator, crossed circles cor- 
respond to matrix elements of the tunneling Hamiltonian, and 
the (red) solid triangles are Cooperon vertices accounting for 
impurity averaging. 

Being interested in low-transparency junctions and re- 
stricting our consideration to the first order in Gi(2), one 
can see that in the case the second electrode is not sub- 
ject to superconducting fluctuations - e.g., is a normal 
STM tip - the only diagram which contributes to the 
tunnel-current is that presented in Fig. [lj. This diagram 
describes the suppression of the tunnel-current due to the 
mechanism of fluctuation rcnormalization of the quasi- 
particlc density of states, discussed above. 

In the absence of magnetic fields, the correlation func- 
tion, Eq.|-(|2|), was already studied in momentum repre- 
sentationU The generalization to the case of a perpen- 
dicular magnetic field can be made by going over from 
the momentum to Landau representation with an appro- 
priate quantization of the Cooper pair motion (see, for 
example, Refs. |j]-|llj). Formally, this corresponds to a 
replacement of the energy associated with the motion of 
the center of mass of a free Cooper pair with momentum 
q by the cigen-encrgy of the Landau state of level m: 
2?q 2 — > lu c (m + 1/2). Here V is the electron diffusion 



coefficient and lo c = AeDH is the cyclotron frequency 
corresponding to the rotation of the center of mass of a 
Cooper pair in a magnetic field H . The integration over 
the two-dimensional momentum in correlator (^) is re- 
placed by a summation over Landau levels according to 
the rule: 



J K 1 m=0 



1 

w c (m+ - 



where M = (T c ot)~ 1 is a cut-off parameter related to 
the elastic electron scattering time r (see Ref. Q for 
details). This transformation is applied to the general 
expression for the correlation function K{l) v ), and one 
findsQ: 



K{iu y ) = K^ {cu y ) + K^ («„) 



(3) 
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[£' m {k + 2 v )-£' m {k)] 
£ m {\k\) 



with a n = e 2 v n T>, Rn being the tunneling resistance of 
the junction and S is its surface area. The function 



£ m (x) = Int + ijj 



1 



Ah 



(4) 



represents the denominator of the fluctuation propagator 
(wavy line in Fig. |l|): 



y n£ m ( x ) 



(5) 



written in Landau representation and describing the fluc- 
tuation pairing of electrons in the normal phase of a 
superconductor over a wide range of temperatures and 
fieldsE. Here t = T/T c0 and h = it 2 /(8-f E )H/H c2 {0) are 
dimensionless temperature and magnetic field normalized 
by critical temperature and the value of second critical 
field respectively, je = 1-78 is the exponential Euler con- 
stant. The cyclotron frequency of a Cooper pair rotation 
in this parametrization is ui c = (16hT c o / tt) . We clarify 
that £' m (x) denotes derivative of the function £ m (x) with 
respect to its argument a;, explicitly given by 



i 



Ah 



(0) 



The two terms in Eq. (|^) correspond to two fluctu- 
ation contributions to the tunnel-current with different 
analytical properties. Below we demonstrate how these 
contributions give rise to the pseudogap maxima and the 
ZBA singularity in the tunneling conductivity in two- 
dimensional disordered superconductors. 



III. COMPLETE EXPRESSION FOR THE 
FLUCTUATION TUNNEL-CURRENT 

We start our analysis with the first term of Eq. (||). 
Since the external frequency w„ enters the expression for 
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FIG. 2. Closed integration contour C in the plane of complex 
frequencies. 



if( reg ) (w„) only via the argument of the analytical func- 
tion £' m (k + 2v) [see Eq. (^)], one can easily perform its 
analytical continuation just by substitution uj v — > —ieV. 
Using Eq. (||) , one finds for the general expression of the 
corresponding current j( reg ) (V): 
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(7) 



The second contribution to the tunneling current is 
determined by 



upper complex plane, where the latter exclude the points 
and iv (see Fig. ||) : 

9 m (w„) = — i^coth (ttz) f m (—iz,u)„)dz . 
c 

The integrals over the vertical line segments become zero, 
the integral over the semi-circle at z = iv is zero since 
fm{k = v,u> v ) = 0, the integral over the semi-circle at 
z = reduces to the residual of coth(7rz). Inverting 
the direction of integration over the line segment with 
Imz = v and then shifting the integration variable as 
z + iuj u /2-kT — > Z\ in the corresponding integral, one 
finds: 

e m K) = + i coth (ttz) 

x [/ m (-iz,w^)-/ ln (-^-^/27rT,w ly )] dz, (10) 

where the "dashed" integral symbol means that the in- 
tegral is performed in the sense of a Cauchy principal 
value. Eq. (|l^) is already an analytical function of cj v 
and one can perform its continuation just by the stan- 
dard substitution u) v — > — iu>. Shifting the variable in the 
second integral again as z — uj/2ttT — > Z2 and using the 
identity 



coth a — coth b = — 



sinh (a — b) 
sinh a sinh b 



one finally finds 



(-icj)-- 



/ m (0, -iw) 



(11) 



with 



if (an) (w„) 
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^^f m (k,u u ), (8) 



m=0 fc=l 



[g;„(2^-fc)-c(fc)] 

&m (k) 



(9) 



Here the analytical contribution is more complex than 
in the case of K ( reg ) (y) , since the frequency w„ is not 
only present in the argument of function (0) but also in 
the upper limit of the sum over k in Eq. (g|). Note, that 
this summation limit can be reduced from v to v — 1 
since f m (k — v,u> v ) = 0. The analytical continuation of 
a function of the form 



E 

fc=i 



f m (k,UJ v ) 



onto the upper half-plane of complex frequencies was per- 
formed in Ref. [@ [see also Ref. [§, equatimi (7.90)]. 
By means of the Eliashberg transformationllj the cor- 
responding sum can be presented as a counterclockwise 
integral over a closed contour C consisting of two hori- 
zontal lines, two vertical lines, and two semicircles in the 



.sinh (w/2T) 



f m (-iz, —iu)dz 



sinh (ttz) sinh7r (z + uj/2t:T) 



Substituting the explicit expression for function 
f m (—iz, —iuj) from Eq. (||) into Eq. ( pd| ) results in 



K^ R (-iw) 



[C(-^)-C(o)] 



ir 3 a n RN „ 

m— 



£m(0) 



(12) 



-2 sinh ( 



[£' m [iz-^)-E' m {-iz)] dz 



\2TJ J £ m (—iz) sinh (nz) sinh n (z + ) 



Eqs. (|2j) and (|12j) determine the second fluctuation 
contribution to the tunneling current 

/(an) (y) 

Let us note that the first term of K^ an ^ R is nothing 
but half of the first summand (with k = 0) of the sum in 
Eq. ([?]) with opposite sign. Technically it would be easy 
to incorporate the latter into K^ reE ' R . However, such a 
procedure would be physically misleading: we will see be- 
low that this k and z independent term in Eq. (^2j) cancels 
the corresponding linear contribution stemming from the 
integral term at small voltages. As a result, the current 
/( an ) (V), determined by the imaginary part of Eqs. (p"2]), 
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does not contain a linear contribution if expanded in pow- 
ers of voltage. This means that it does not contribute 
to the magnitude of the differential tunnel conductiv- 
ity at zero voltage a^l (T, H, V = 0) = dI^/dV\ v=0 , 
which is the easiest quantity to measure in experiments. 
Nevertheless, it contributes to the current- voltage char- 
acteristics at finite voltages and, as we will see below, 



can noticeably manifest itself even at very low voltages 
eV ~ T - T c0 as a ZBA. 

Adding /( r °s) (V) and /< an ) (V) one finds the general 
expression for the fluctuation contribution to the tunnel- 
current, which is valid in the complete phase diagram 
beyond the H c2 (T) line: 



I (fl) (t, h, V) = 7 (rcg) + I (an) 
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dz 



m=0k=0 v ' m=0 

Re£ m (iz) [Re£^ (iz - ieV) - Re£^ (iz)} + Im£ m (iz) [Img^ (iz - ieV) + lm.£' m (iz)} 
sinh (irz) sinh [tt (z - eV/2irT)} [Re 2 £ m {iz) + Im 2 £ m (iz)] 
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FIG. 3. (Color online) Illustration of the t-h phase diagram 
with the parameter space location of all plots in this paper. 
In each plot the parameter points (A; or Bf), parameter line 
cuts (Li), or line groups (LLi) are given. 

Eq. ( |l3|) is the main result of this work. The first term 
/0 c g) (y) has been studied in detail for diffcrejrtdLmiting 
cases using different approaches: close tc^T c oJj'Qoll3: (i) 
in a wide temperature range in zcrOjfipldEI, or (ii) close to 
T c q in magnetic fields H <§; H c2 (Q)}l3. The current con- 
tribution /( an ) (V) has been omitted in all these works 
based on the "standard" argument that the zero fre- 
quency bosonic mode (which traverses through the prop- 
agator) is singular in the vicinity of the transition. How- 
ever, it is known that this argument sometimes works 
(e.g., in the .case of the Maki- Thompson contribution to 
conductivity!!!!) , but also sometimes fails (e.g.,r-fcr the 
Aslamazov-Larkin contribution to conductivity^) . In 
our case this argument turns to out be correct only for 
very small voltages. The reason being that voltage itself, 
together with temperature deviations from the transition 
point and finite magnetic fields, drives the system away 
from the immediate vicinity of the transition, which in- 
validates the argument regarding the dominance of the 



zero frequency bosonic mode. 

We present several plots of the tunnel-current and the 
tunnel conductance. Since they depend on three parame- 
ters: t, h, and v, only lines or planes in the full parameter 
space are presented as line or surface plots. We included 
Fig. ||, showing all parameter points and lines in the t-h 
phase diagram for all following figures. The critical field 
line, h c2 (t), separating the superconducting (SC) and the 
normal fluctuation region (SF) is defined by £q(0) = 0. 
Each figure caption refers to these parameter locations. 
Fig. shows the behavior of I^ & > (t, h, V) near T c q and 
H c2 (0). 

In the following we will carefully analyze the effect 
of superconducting fluctuations in the whole phase di- 
agram. We start our discussion with the regular contri- 
bution /( ro s) (V) and then elucidate the important role 
of the anomalous contribution /^ an ) (V"), which was ne- 
glected in literature so far. 



IV. ANALYSIS OF THE ASYMPTOTIC 
BEHAVIOR OF J (rcg) (V) 

Close to T c q and for sufficiently weak magnetic fields 
H <C H c2 (0), the most singular term in Eq. (f?j) arises 
from the zero frequency bosonic mode k = 0, when the 
propagator has a pole at e = and 

£ m (0) = e + 2h (m + (14) 

with e = lnt w t - l^las reduced temperature. The 
summation over Landau levels can be performed in terms 
of polygamma-functions, tp( n '(x), and one finds an ex- 
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FIG. 4. (Color online) Total tunneling current close to T c q (left) and near h C 2{0) (right) at various temperatures depending on 
the dimensionless voltage v = 2eV/ABCS- The insets show the regular and anomalous contributions at the respective lowest 
temperature separately. As one can see, the anomalousnart has a nonlinear component near v = 0. The current is normalized 
to Jo = eT c o5'/(o"„J?jv). (left) parameter points in Fig. [j are A1-A4, inset Ai, (right) parameter points are B1-B3, inset B\. 



pression valid for any combination of e and d«l 
eTS 
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(15) 



Eq. ( |l5|) reproduces the results of Refs. [0 and The 
corresponding contribution to the tunneling conductance 
is 
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In the region of high temperatures T 3> T c o and zero 
magnetic field we restrict our analytical consideration to 
the fluctuation contribution to the differential conductiv- 
ity at zero voltage. Performing an integration instead of 
a summation in Eq. ^ one finds 



a^(0,t > 1) 



Se 2 



Ait 2 o n R 



In- 



to 



T c0 r 



N 



lni 



which is again in complete agreement with Rcf. 0. 

Close to the line H C 2 (t) and for sufficiently low temper- 
atures t <C h C 2 (t) the lowest Landau level approximation 



(LLL) holds. The corresponding propagator (with quan- 
tum number m = 0) has a pole structure and Eq. (Q) 
acquires the form: 



So (k) = h 



ir 2 tk 
4h~ 



c:2 



(17) 



with h(t) = (H - H c2 (f)) /H C 2 (t). Keeping only the 
m = term in Eq. (Q), one can write 



1^ [V,t<h c2 (t)] 



2eT c oSh 
■K 3 a n R N 



T 1 ir 2 tk 
k=Q n 4/i c2 (t) 



. (18) 

The imaginary part Im£g (fc — ieV/irT) can be explicitly 
written using Eq. (^) in the limit t -C h C 2(t) and the 
asymptotic behavior of ip' (\x\ 1) ~ 1/x: 



Im £' n [ k — ■ I = 

ttT ) 2ttT 
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The summation in Eq. (|l8]) can then be performed ex- 
actly in terms of polygamma- functions, i.e., using 



E 

fe=0 



J 1 _ 1 (o) - %b ((3 + ij) 

: + ot (k + f3) 2 +j 2 ~ 7 fB + ij-a 



which gives an expression for the regular part of the fluc- 
tuation current valid for low enough temperatures along 
the line h C 2(t): 



1^ [V t ,t<&h c 2 (t)]=- 



2eST c0 h v t 
ir 3 a n RN 1 + v 2 



Here, we introduced the dimensionless voltage 
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(20) 
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which defines the characteristic scale of the con- 
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sidcrcd domain of the phase diagram. We stress, that this 
scale depends on temperature via the parameter ft c2 (t). 
Close to H C 2 (0), in the region of very low temperatures 



t -C ft, the argument of the ^-function in Eq. (20) be- 
comes large despite the smallness of ft, and the ^-function 
can therefore be approximated by its asymptotic expres- 
sion. One gets 



/(rcg) ^ 
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(21) 



with Abcs = kT c q/^e being the value of BCS gap. The 
characteristic scale where the maximum of the tunnel 
conductance appears at these low temperatures is v = 
2eV/A BCS ~ 1, i-e. 



eVmax ~ Abcs- 



(22) 



In the region of high fields H 3> H c i and low temper- 
atures, the asymptotic behavior of the tunneling current 
can be studied in complete analogy to the case of high 
temperatures and weak fields. The sums in Eq. (0) can 
be approximated by integrals, which gives for the value 
of the differential conductivity at zero voltage: 



cr (rog) (0,ft>l) = - 
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One can see that this dependence is exactly the same as 
that one in the case of high temperatures with reversed 
roles of the reduced temperature and the reduced held. 



V. LOW VOLTAGE BEHAVIOR OF J (an) (V) 

In the low-voltage limit, V — > 0, the general expression 
can be expanded in small eV. We start 
with the first order term of that expansion, where one can 
assume V = in the argument of integrand function and 
obtain 
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In the region of temperatures close to the transition 
temperature T c o and along the transition line for temper- 
atures t <?C ft C 2 (tY the propagator has a simple pole struc- 
ture [see Eqs. (jlj) and (0)] and the integral in Eq. ( |23] ) 
can be calculated analytically. Performing this integra- 
tion one finds that the second term of Eq. (]2^) exactly 
annihilates the linear part of the first term. This fact 
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FIG. 5. (Color online) Regular and anomalous contributions 
to the tunneling conductance close to T c o (top) and at low 
temperatures near h C 2(0) (bottom). The regular part is pre- 
sented by a solid line (red), the anomalous by crossed line 
(green) and the their sum, i.e., the total fluctuation contribu- 
tion, is shown by a dashed line (blue), (top) parameter point 
in Fig. | is A\, (bottom) B\. 



justifies the static approximation (zero frequency bosonic 
mode) made in Rcfs. [z| [H] . Yet, this static approx- 
imation turns out to be valid only for very low voltages. 
Expanding the integrand in Eq. dl3| ) to higher orders in 
voltage reveals an unexpected result. One can see that 
the voltage V enters the integrand of Eq. ( |l3|) in two 
different places: in the argument of £ m (iz — ieV) in the 
numerator and in the argument of sinh (~kz — eV/2T) in 
the denominator. The expansion of £ m (iz — ieV) results 
in the appearance of a weakly voltage-dependent term 
of the order of O(V 3 /T^ ) in /( an ), while, as one can 
easily verify, the expansion of sinh -1 (ttz — eV/2T) after 
integration leads to a very singular correction 
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At zero magnetic field, the dip in the corresponding 
tunnel conductivity rapidly develops on the scale eV ~ 



7 




T — T c q, while at zero temperature close to H C 2 (0) this 
happens on the scale eV ~ Abcs' 1 - 
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FIG. 7. (Color online) Temperature and voltage dependence 
of the tunneling conductance due to superconducting fluctu- 
ations above T c o in zero magnetic field. The corresponding 
parameter line in Fig. ^|is L\. 
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FIG. 6. (Color online) Smearing of the low-voltage singular- 
ity of the tunneling conductance with increasing temperature 
close to T c o (top) and near /i C 2(0) (bottom), (top) parame- 
ter points in Fig. ^ are A\-A±, (right) parameter points are 
Bx-Ba. 

The effect of both fluctuation contributions, /( reg ) (V) 
and /( an ) (V), on the tunneling conductance is demon- 
strated in Fig. |^. Similar behavior can be observed along 
the whole line H C 2 (T) . The singularity in the low volt- 
age behavior of tunneling conductance rapidly smears out 
when moving away from the transition line or increasing 
the temperature (see Fig. ||). 

VI. NUMERICAL ANALYSIS 

The temperature, magnetic field, and voltage depen- 
dencies of the tunneling conductance due to supercon- 
ducting fluctuations, calculated numerically based on 
Eq. (|i3|), are presented in Figs. 0- 1| as surface plots. The 
numerical procedure to calculate the fc-sum of the first 
term needs to take into account its relatively slow conver- 
gence. Therefore it is calculated explicitly up to a thresh- 



FIG. 8. (Color online) Magnetic field and voltage dependence 
of the tunneling conductance due to superconducting fluctu- 
ations at t = 0.05 above h C 2(0) = 0.69. The corresponding 
parameter line in Fig. |§] is L2 ■ 

old at which the sum can be replaced by an integral and 
the polygamma functions by their asymptotic behavior, 
(here we use as thrcshold-fc, the value &m at which the 
argument of the function £ m reaches 1000). The "rest"- 
intcgrals are calculated with inverse integration variable 
using a Gauss-Lcgcndre method. The second term re- 
quires a careful treatment of the two integrable poles, 
which is done by analytical calculation of the residuals 
in a small interval around them, where the denomina- 
tor is linearized. Also the numerical integration outside 
the pole intervals is done by using adaptive integration 
point distances. The overall behavior of both terms of 
the tunnel-current results in a pronounced pseudo-gap 
structure of the conductance near the superconducting 
region. It is the non-linear anomalous term of the tunnel- 
current which is responsible for the fine structure (local 
maximum) at the center of the gap, the ZBA. 

At this point it is worth mentioning that another sharp 
fine structure of tunnel conductance which should oc- 
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cur in the same scale eV ~ T — T c q was predicted in 
Rcf. Q. This structure appears due to interaction of 
fluctuations as the second order correction in Ginzburg- 
Levanyuk number Gh 2 ) (but still in first order in the bar- 
rier transparency). This contribution has an interference 
nature (analogously to Maki-Thompson process) and, in 
contrast to the discussed above nonlinear contribution 
CT (an) ( e y <<; T _ Tcq) „ Gi (2) [eV/(T - T c0 )] 2 , diverges 

at zero voltage as Gi 2 2) [T c0 /(T - T c0 )] 2 ln[(T - T c0 )/eV}. 
Such divergency, in complete analogy to Maki-Thompson 
cont ribjut ion, is cut off by any phase-breaking mecha- 
nismLM 

Analyzing the surface plot representation of the ex- 
perimental results of Ref. 0, obtained at temperature 
close to T c o, one notices their striking similarity to the 
theoretical surfaces presented in Fig. [?]. Indeed, the au- 
thors of Ref. H mentioned the agreement of their results 
with the theoretical prediction of Ref. 0. Fig. || shows 
how the corresponding surface transforms at low temper- 
atures and strong magnetic fields close to H c2 (0). 
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FIG. 10. (Color online) Zero-bias tunnel conductance as 
function of magnetic field near the critical temperature for 
t — 1.0, 1.025, 1.05, 1.1, 1.5. These correspond to the parame- 
ter line group LLi in Fig. 0. 
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FIG. 9. (Color online) Zero-bias tunnel conductance, 
a(t,h,v — 0), as function of temperature at fields 
slightly above the critical zero-temperature field for h = 
0.695,0.7,0.705,0.71,0.75. These correspond to the param- 
eter line group LL2 in Fig. |^. 

It is interesting to note that the behavior of the gen- 
eral expression ([13]) clearly shows growth of the fluctua- 
tion effects in the domain of intermediate temperatures 
and magnetic fields, beyond the immediate vicinity of 
T c o and iJ C 2(0), see plots of the zero-bias tunnel conduc- 
tance a(t,h,0) = cr( rog) (t, h, 0) in Figs. | and [[(]. In 
Fig. [y] one can see the evolution of the pscudogap near 
the h C 2(t) line (slightly offset by a factor 1.1, see cap- 
tion), exhibiting a deeper suppression for intermediate 
temperatures and fields. This fact is in agreement with 
the general ideas of the theory of fluctuations establish- 
ing the growth of fluctuations strength (characterized by 
the Ginzburg-Levanyuk number) as one moves away from 
the extreme points [T c q and H C 2(0)] of the curve H C 2(T) 



(see chapter 2 of Ref. ||). 



0.0 




FIG. 11. (Color online) Voltage dependence of the tunneling 
conductance due to superconducting fluctuations along the 
h C 2(t) plotted for (t, h) — s(x,h C 2(x)) with x £ [0; 1] and 
distance s = 1.1. The corresponding parameter line in Fig. ^ 
is L3. 



VII. DISCUSSION AND COMPARISON WITH 
EXPERIMENTS 

We have calculated the effect of SF on the tunnel- 
current characteristics for a SIN junction of a two- 
dimensional disordered superconductor and a normal 
metal electrode. In contrast to the common believe that 
the shape of dtun, Eq. (|IJ), presents an almost direct im- 
age of the energy deoendence of the density of states of 
a measured samplccJ, one notices, that this is not the 
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FIG. 12. (Color online) Theoretical curve of the fluctuation 
correction to the single particle DOS, dv^) i > versus ener gyi E, 
for 2D superconductors above critical temperature for tem- 
peratures t = 1.05, 1.2, 1.5, 2.0. The inset shows more details 
for temperatures close to T c near the maximum. The corre- 
sponding parameter points in Fig. ^ are A1-A3. 

case for the SIN system studied here. Indeed, the ef- 
fect of fluctuation Cooper pairing above T c results in the 
depletion of the electron density of states at the Fermi 
level [E = 0) and its growth at the characteristic ener- 
gies of the Cooper pairs, E = T — T c . This phenomenon 
was considered theoretically long ago (sec Rcf. |2 andH]) 
and the corresponding curves for b~v)$l (E) in the absence 
of magnetic fields and not too far away from the Fermi 
level are presented in Fig. n3. Comparing these curves 
to our cr^jj (V, T) (see Figs?~|| and ^) reveals an appar- 
ent discrepancy, which is a result of the dramatic energy 
dependence of Sv^J (E) at small energies (E <T). 

Note, that the statement concerning the similarity of 
the differential tunnel conductance versus voltage and 
the energy dependence of the density of states, is based 
on the common assumption that at low temperatures 
the latter changes on the scale E 3> T and drip /BE ~ 
cosh -2 E/2T acts as a delta- function in Eq. (|l|). This is 
exactly the case of the ZBA occurring in disordered met- 
alsEJ, where the correction to the density of states due to 
inter-electron interaction varies on the scale E ~ t _1 3> 
T and 5<r[^ A) (V) ~ cM Coulomb ) (V). 

However, the situation is just the opposite for the 
single-electron density of states which is renormalized 
by superconducting fluctuations. Here the function 
<5i^ 2 -j (E, T) strongly varies and even changes its sign on 
the scale E ~ T — T c o -C T, while vanishing at large ener- 
gies. In the energy interval, E £ [0; 0(T)], the derivative 
driF/dE remains practically equal to one. Extending 
the approximation dnp/dE = 1 to the complete inte- 
gration domain E £ [0;oo), one would obtain zero due 
to the sum rule Sv^) (E, T) dE = (the latter re- 
flects the conservation of the number of states). In or- 
der to obtain a finite conductivity, one needs to ana- 



lyze Eq. (pi) more carefully, taking into account the high- 
energy tail of (E,T), which is beyond the approx- 
imations used in Refs. |^ and |3j. Following the method 

proposed in Ref. 0), we calculated a[^ D (V, T) diagram- 
matically, without ever referring directly to the density 
of states, which consequently avoids the above mentioned 
limitations. As a result of our detailed analysis of the 
tunnel-current, the wide bell-shaped form of the tunnel 
conductance with peripheral maximum at eVmax ~ ksT C Q 
and sharp peak occurring at zero voltage, the ZBA, was 
found (see Figs. || and ||). 

In-jjpC|Cnt experimental works on Bi-2212 and Nb sam- 
plescBEil, very interesting and carefully measured results 
on the temperature and magnetic field dependence of the 
differential tunnel conductivity arc presented. Unfortu- 
nately, a direct comparison of our results with those ex- 
perimcftts-is impossible at this point, since the authors 
of Refs.EK3 investigated SIS junctions, while we consider 
SIN (or STM) systems here. Nevertheless, a careful study 
of Refs. [22] and reveals several similarities between 
the experimental finding and our prediction. First, the 
authors always observe a peak in the fluctuation region 
(T > T C (H)), which corresponds to energies smaller than 
2A BCS {T = 0) (Figs. 11 & 12 of Ref. g]). Secondly, 
very characteristic is the presence of a "crossing point" 
in the temperature dependencies of the differential tun- 
nel conductivity (see Figs. 4 & 5b of Ref. |2^]), analogous 
to that one of our Fig. Finally, the experiments also 
showed a very strong dependence of the tunnel magneto- 
conductivity on voltage close to V = (see Figures 6 
d.c and 7 of Ref. p^|), which could be identified as the 
predicted nonlinear anomalous contribution [the second 
term of Eq. ([j"3|), which is rapidly smeared out when one 
moves away from the transition line or increases temper- 
ature (see Fig. ||)]. However, in view of the SIS-nature 
of the junctions used in the experiment, the observed 
ZBA could also be attributed to presence of a residual 
Joscphson current. To the best of our knowledge this 
is the only case where the existence of a ZBA in the 
tunnel-conductivity close to the critical temperature is 
reported. A definitive identification of its origin requires 
further studies. 

The difficulty in observing the predicted ZBA due to 
can be related to the strong dependence of this 
ZBA on inhomogeneities in the barrier and local varia- 
tions of the critical temperature. Regrading the above 
mentioned second order correction due to the interaction 
of fluctuations, it was never reported to be observed ex- 
perimentally. 
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